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Describing matter at near absolute zero temperature 
requires understanding a system's quantum ground state 
and the low energy excitations around it, the quasiparti- 
cles, which are thermally populated by the system's con- 
tact to a heat bath. However, this paradigm breaks down 
if thermal equilibration is obstructed. While such non- 
equilibrium behavior may be expected in disordered and 
frustrated quantum systems (like for instance quantum 
spin glasses long-range Josephson junction arrays in 
a frustrating magnetic field 2], or self-generated mean- 
field glasses 3]), it is non-obvious that it may exist in 
clean systems with only local interactions and without a 
complicated distribution of energy levels. In this paper 
I present solvable examples, three dimensional general- 
izations of Hamiltonians proposed for topological quan- 
tum computing, that have solely local interactions, no 
quenched disorder, and relaxation rates akin to those of 
strong structural glasses. Therefore, in these systems the 
topologically ordered ground states are not reached when 
the temperature is reduced to absolute zero. 

Topological order and quantum number fractionaliza- 
tion are some of the most remarkable properties of sys- 
tems of strongly interacting particles. Some phases of 
matter, in contrast to common examples like crystals and 
magnets, are not characterized by a local order parameter 
and broken symmetries. Instead, as shown by Wen 
some quantum phases are characterized by their topo- 
logical order, such as the degeneracy of the ground state 
when the system is defined on a torus or other surface 
of higher genus. These topological degeneracies cannot 
be lifted by any local perturbation. Topological order 
and quantum number fractionalization are intimately re- 
lated, and much effort has recently been directed at these 
exotic properties, for they may play a role in the mech- 
anism for high-temperature superconductivity H 0, @ • 
Also, the robustness of a topological degeneracy to local 
noise due to an environment is at the core of the idea 
behind topological quantum computation, as proposed 
by Kitaev Interestingly enough, strong correlations 
that can lead to these exotic quantum spectral proper- 
ties can in some instances also impose kinetic constraints, 
similar to those studied in the context of classical glass 
formers [3 El 13 111 111 13 13 E3 
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FIG. 1: (a) Cubic cell of an fee lattice. The centers of the 
six faces form an octahedron, with its sites labeled from 1 
(topmost) to 6. In addition to the set of octahedra formed by 
the face centered sites, there are three more sets of octahedra 
that can be assembled from sites both on faces and on cor- 
ners of the cubic cells, totaling 4 such sets. Six-spin operators 
are defined on these octahedra using the <r x ' y ' z components of 
spin on each vertex as described in the text, 
(b) Centers of 6 octahedra cells that share a spin, which re- 
sides at the site / shown at the center. The x,y, or z labels 
sitting at the centers of the octahedra show which spin op- 
erator <Tj' y ' z flip their Oi eigenvalue. Acting with any of the 
operators <r^' y ' z always flip the eigenvalues Oi of exactly four 
octahedra. 



The possibility of glassiness in pure strongly corre- 
lated quantum systems with solely local interactions is 
demonstrated by studying the following exactly solv- 
able example. A model displaying strong like glassi- 
ness is constructed on a three-dimensional (3D) face- 
centered cubic (fee) Bravais lattice, spanned by the prim- 

itive vectors a x = (-*=,-^,o), a 2 = (fj,-^,-^), and 

a 3 = ^-^5, 0, -j=j . Each site can be indexed by i, j, k G Z, 

denoted by a superindex / = (i,j,k). At every lattice 
site / one defines quantum spin S — 1/2 operators a*, 
<7j, and erf. 

The fee lattice can house sets of octahedra: the simplest 
one to visualize is the one assembled from the centers of 
the six faces of a cubic cell, and is shown in Fig. Ofa). 
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In addition to this simple set, there are three more sets 
of octahedra that can be assembled from sites both on 
faces and on corners of the cubic cells, totaling 4 such 
sets, which we label by A, B, C and D. 

It is simple to see that the total number of octahe- 
dra equals the number of spins: each lattice site I is 
the topmost vertex of a single octahedron. Define Pj as 
the set of six lattice points forming the octahedron with 
site I at its top. The six vertices are indexed by J n (I), 
for n — 1,...,6, with one of the vertices J\{I) = I- 
The six labels are assigned in such as way that the pairs 
{ Ji, J a], { J2, J5}, { J3, Jq\ are diagonally opposite sites 
from one another, and this number labeling is illustrated 
for a single octahedron in Fig. Qfa). From the one-to- 
one relation between a site / and an octahedron Pi, we 
can also partition the lattice sites into four sublattices 
A, B, C and D. 

Now define the operators Oj as 

°i = ct j 3 (/) CT J 4 (/) a Mi) ■ (!) 

This construction generalizes to 3D the plaquette inter- 
actions defined for planar 2D lattices by Kitaev 9] and 
Wen [18| . These operators commute, [Oi,Op] = for 
all pairs /, /'. It is simple to see how: two octahedra Pjr 
and Pp can either share 0,1, or, at most, 2 spins. If they 
share spins, they trivially commute. If they share 1 
spin, the component (x, y or z) of er for that shared spin 
coincides for both Oj and Op (the two octahedra touch 
along one of their diagonals). If they share 2 spins, the 
components of <x used in the definition of Oi and Op 
are different for both spins, there is a minus sign from 
commuting the spin operators from each of the shared 
spins, and the two minus signs cancel each other. 
Consider the Hamiltonian 

& = ~ y £o x . (2) 

Because the Oi all commute, the eigenvalues of the 
Hamiltonian can be labeled by the list of eigenvalues 
{Oi} of all the 0j. Notice that Oj = 1, and so each 
Oi = ±1. In particular, the ground state corresponds to 
Oi = 1 for all I. 

Because the number of spins equals the number N of 
sites and of octahedra, one may naively expect that the 
list {0/ = ±1} exhausts the 2 N states in the Hilbert 
space. However, there are constraints that the Oi satisfy 
when the system is subject to periodic boundary condi- 
tions (compactified). One can show that 

Y[oi=noi=noi=i[oi=t. (3) 

IeA l£B JeC IGD 

There are four constraints; therefore there are only 2 
independent {Oi = ±1}. This implies, in particular, 
that there is a ground state degeneracy of 2 4 = 16. The 
ground state degeneracy is not associated with a symme- 
try; in particular, it is easy to show that (<Jj' y ' z ) = 0. 



This is a topological degeneracy, and the eigenvalues 
T a = ±1 (a = 1,2,3,4) of a set of four non-local (topo- 
logical) operators T a are needed to distinguish between 
the 16 degenerate ground states. 

The operators T a can be constructed as follows. Let 
'Pi = {I\j + k = 1} be a set of points along a horizontal 
plane. Notice that each such plane contains sites in only 
two of the four sublattices A, B, C, D. For example V\ C 
A U B and V 2 C C U D. Define Ti = UiePinA a h T 2 = 

ri/e-PinB /' T3 — Uiev 2 nc' J h ^4 = Hiev 2 nD °/ ■ ^ 
is simple to check that [T a , Op) = for all a and I, and 
the T a trivially commute among themselves. Hence the 
four eigenvalues ^,2,3,4 = ±1 of ^2,3,4 can distinguish 
the 16 degenerate ground states. 

The spectrum of the Hamiltonian Eq. (J2J) is that of a 
trivial set of N — 4 free spins, determined by the list of 
eigenvalues {Oj = ±1} of all the Oi, subject to the con- 
dition Eq. ©: E{ 0l } = — jE/^' Excitations above 
the ground state (Oi = 1 for all I) are "defects" where 
Oi = — 1 at certain sites I. The equilibrium partition 
function is given by Z = 16 J2{Oi=±i} e 5/3/l ^ J °' ■ At 
thermal equilibrium at temperature T, the thermal aver- 
age (Oi) = tanh and the concentration or density of 
Oi = —1 defects is c = i (l — tanh ^). Notice that we 
have encountered an analogous situation to that in the 
classical spin facilitated models 01 > m particular the pla- 
quette models displaying glassy dynamics p^ll^lla[l7l| : 
the thermodynamics is trivial in terms of non-interacting 
defect variables. To study the approach to such an equi- 
librium state, the coupling of the system to a bath of 
quantum oscillators must be introduced. 

The Hamiltonian of the system plus bath of oscillators 
can be formulated as |Hj, |2(j 

7~i = H + iibath + -^spin/bath 

where H is defined in Eq. (J2J, the bath -ffbath contains 
a family of harmonic oscillators a\j,a\ l for each site, 
and 

l,a A 

where the g a are the generic coupling constants for each 
of the three components (a = 1,2,3 or x, y, z) of the 
spins. 

Although the spectrum of H in Eq. J3J) is the same 
as that of free spins in a uniform magnetic field h, the 
variables 0/ for different octahedra P x cannot be in- 
dependently changed, as opposed to spin variables in 
a free spin model in a field h. The bath couples to 
the physical degrees of freedom, the spins crj. Acting 
on a site I E Pi with one of the operators (Jj,er~, or 
cr| flips or not the eigenvalue Oj depending on whether 
oy y ' z Oi = tOi o-j' y ' z , respectively. However, the spin 
(Tj is shared by six neighboring octahedra, and thus one 
cannot change the eigenvalue of Oi without changing the 
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eigenvalues Op of some of the neighbors by the action of 
the local spin operator. 

If integrated out, the bath degrees of freedom lead to 
a non-local in time action and to dissipation effects. In- 
stead of working with the dissipative action, let us fol- 
low the time evolution of the system plus bath, and look 
at the possible evolution pathways of the quantum me- 
chanical amplitudes of the system plus bath degrees of 
freedom. (Yet another alternative is to work within the 
von Neumann density matrix formalism [2l| . and follow 
the time evolution of the matrix elements). After evolu- 
tion for time t from some initial state, the system is in a 
quantum mechanical superposition 

!*(*)> = E r {T..o,}(t) KT«,Oj}> ® |T {Ta ,o,}W) , 

{T a ,0 I= ±l} 

(5) 

where \^{T a ,Oi}(t)) is some state in the bath Hilbert 
space with norm one. The fact that the bath degrees of 
freedom couple to single quantum spins 07 [as in Eq. 
enters the problem through the permitted channels for 
dynamically transferring amplitudes among the ^{T a ,Oi}- 

The processes that redistribute or transfer amplitude 
among the Tpp Oj} correspond to different orders in per- 
turbation theory on the g a system-bath coupling. There 
is also a thermal probability factor coming from the bath 
that depends on the difference between the initial and 
final energy -E{Oj} — — f J2i Oi 01 * ne system. The sim- 
plest class of paths is a sequential passage over states 
connected through first order in g a processes; this is a 
"semi-classical" type trajectory. 

Through the action of a local <7~' y ' z operator, exactly 

4 of the 6 octahedra operators Oj sharing spin 7 are 
flipped. The reason is that the 6 octahedra operators can 
be divided into 3 groups of 2 octahedra having in their 
definitions Eq. respectively, the x,y, and z component 
of spin operator at the shared site. This is illustrated in 
Fig. IHb). Acting with either of the three components 
of the spin operator cr~' y ' z on this shared site will flip 
exactly 4 out of 6 defect variables Oj. Hence, a sin- 
gle defect cannot be annihilated in this process. Defects 
disappear from the system only through recombination. 
This multi-defect type dynamics makes it difficult for the 
system to relax to equilibrium, exactly as in the kinet- 
ically constrained classical models [2 0, 0, . For 
example, if the temperature is lowered, in order to de- 
crease the defect density, either four defects must come 
together and annihilate (4^0 decay), or three defects 
become one (3 — > 1 decay). Moreover, single defects 
cannot simply diffuse through the system; that would 
require flipping only two neighboring octahedra, but in- 
stead four are always flipped. To move, an isolated single 
defect must first decay into three defects (1 — ► 3 produc- 
tion) because of the multi-defect dynamics, then a pair 
can diffuse freely (2 — > 2) and recombine with another 
defect through a 3 — ► 1 decay process. Because of the 
initial 1 — > 3 production process, there is an energy bar- 
rier of 2h to be overcome. This activation barrier leads 



to recombination/equilibration times 



scq. 



T exp(2h/T) 



that grow in an Arrhenius fashion as temperature is low- 
ered (to is a microscopic time scale). 

What about quantum tunneling? Amplitude can be 
transfered from some initial to some final state via virtual 
processes, in which the number of defects is larger in the 
intermediate (virtual) steps. Virtual processes of nth or- 
der involve a product of n spin operators, T = Yi s =i a f' ■ 
For a single defect to disperse through quantum tun- 
neling, an J- operator that flips only two octahedra is 
needed. However, one can show that any T will flip at 
least four octahedra (as opposed to Kitaev's and Wen's 
models, in which two defects stand at the endpoints of 
"strings" , here four or more defects lie at the corners of 
"membranes"). Although defects cannot disperse, tun- 
neling still contributes to defect annihilation and to de- 
fect pair motion. In perturbation theory, a process in 
which a defect pair separated by a distance £ can hop by 
a lattice spacing has an amplitude of order (g/hp (no- 
tice the energy denominator h) and defect annihilation 

has an amplitude (g/h)^ . Hence virtual processes are 
suppressed exponetially in £, and if the system were to 
equilibrate at temperature T, where the typical defect 
separation is £ = c -1 / 3 ~ e h / 3T , the characteristic time 
scale 



T exp 



Hh/g) e h ^ T 



would grow extremely quickly as the temperature is low- 
ered. What we learn from this simple estimate is that 
quantum tunneling is less effective than classical sequen- 
tial processes in thermalizing the system. This is coun- 
terintuitive to the notion that at low temperatures quan- 
tum tunneling under energy barriers remains an open 
process while classical mechanisms are suppressed due to 
high thermal activation costs. The reason for the par- 
ticular quantum freezing in this system is simple: as the 
distance between defects increase at lower temperatures, 
the barrier widths increase, which debilitates tunneling. 
In passing, we note that in a finite system of size L, one 
must replace £ by L in the estimation of the recombina- 
tion/equilibration times, t tU n. ~ To exp \ln(h/g) L]\ this 
time scale is also of the order of that for tunneling be- 
tween two topological ground states in a finite system of 
size L 0. 

Because i soq . and t tun . grow rapidly as the tempera- 
ture lowers, the system will fall out of equilibrium at low 
temperatures, and physical correlation functions will not 
be those simply computed in the framework of equilib- 
rium quantum statistical mechanics. The simplest corre- 
lation function that illustrates this point is the one-point 
function related to the time-dependent spatial concentra- 
tion of defects pi(t) = |[1 — (Oi(t))]; let us find how it 
approaches, as a function of time, the asymptotic equi- 
librium value Cf = i ^1 — tanh when the tempera- 
ture is, say, reduced by half from 71 to 71 = 71/2. The 
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mechanism for equilibration is diffusion-annihilation of 
defects. We have argued that defects are always flipped 
in quadruplets, and single defects cannot freely diffuse 
without generating more defects. Defect pairs, however, 
are free to diffuse quickly, so we can reduce the problem 
to an effective reaction-diffusion 22] of the A + A ^ 
type for the single defects facilitated by the pair motion. 
The quantum average over the state Eq. (JSJ can restore 
translational invariance of the density pi(t), so spatially 
homogeneous initial densities remain homogeneous un- 
der time evolution, pi(t) = c(t), hence the dynamics for 
this problem is controlled by the simple rate equation 
c(t) = —k(c 2 — Cf ) , with the kinetic rate coefficient k 
directly proportional to the defect diffusion constant at 
temperature Tf, from which it follows that k oc l/t seq .. 
At long times c(t) — Cf oc Cf exp(— 2cffc t), from which 
we extract the time constant for the relaxation of the 
one-point correlations to be ri pt = t scq ./2cf. Notice that 
the relaxation time r lpt for the one-point correlation is 
longer than t scq . because the annihilation rate is reduced 
for low densities of defects. This enhancement must be 
cut off when the density of defects is of order 1/L 3 , in 
which case ri pt ~ I? i soq .. The relaxation time is just 
polynomial in the system size L, so if t scq . saturated to 
a constant value as temperature is lowered, the system 
would not behave as a glass, ft is the Arrhenius form of 
t scq . that causes the glassy behavior. 

The fact that the system presented above is exactly 
solvable helps to understand the origin of its glassy be- 
havior, but it is not a necessary ingredient. To illus- 
trate this point, one can simply add a perturbation 
SH = J^ia^a a i- For <^a/h less than order unity, 
this interaction can be analyzed in perturbation theory 
similarly to the arguments above for the ratios g a /h. The 
perturbation will give the defects some mobility, but that 
becomes exponentially small as the defects grow apart. 
Indeed, these arguments can be generalized for any local 
perturbation written in terms of the physical spins <xj, 
as long as the coupling constants are small compared to 



the gap h. 

To summarize, the essence of why this quantum sys- 
tem is glassy is the following. The thermodynamics is 
best described by working in the basis of eigenstates or 
defects \{T a ,Oi}); however, upon acting on these states 
with physical spin operators crj' y,z , single defects can nei- 
ther be annihilated nor simply moved around (diffused). 
The lack of defect diffusion in these glassy systems is 
protected by the fact that any physical spin operators 
must flip quadruplets, not pairs, of defect variables. The 
system can only relax by multi-defect real processes that 
are thermally suppressed or else by virtual processes that 
involve quantum tunneling of increasingly large objects 
as the defect density is reduced at low temperatures. 

Many elastic, thermal, electronic, and magnetic prop- 
erties of classical glassy material systems are conse- 
quences of these materials' being out of equilibrium. 
Such properties can be tailored according to preparation 
schemes - for example, by controlling cooling rates. In 
contrast, because of the difficulties in studying real-time 
dynamics of strongly interacting quantum systems cou- 
pled to a thermal bath, very little is currently known 
about properties of quantum matter that can be en- 
gineered by keeping systems out of equilibrium. In a 
broader scope, the main result of this work is that it 
presents a concrete example of a solvable toy model which 
shows without arbitrary or questionable approximations 
that a pure quantum system with only local interactions 
can, indeed, stay out of equilibrium. This result supports 
the possibility that there may be material properties due 
to non-equilibrium glassy behavior in quantum matter. 
It also suggests a new design constraint for topological 
quantum computing: that the ground state degeneracy 
is protected while the system is still able to reach the 
ground states. 
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